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Abstract 



A theoretical study of the specific heat C(T) as a function of temperature 
in Fibonacci magnetic superlattices is presented. We consider quasiper iodic 
structures composed of ferromagnetic films, each described by the Heisenberg 
model, with biquadratic and bilinear coupling between them. We have taken 
the ratios between the biquadratic and bilinear exchange terms according 
to experimental data recently measured for different regions of their regime. 
Although some previous properties of the spin wave specific heat are also 
reproduced here, new features appear in this case, the most important of them 
being an interesting broken-symmetry related to the interlayer biquadratic 
term. 
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I. INTRODUCTION 



The study of the properties of magnetic multilayers has been one of the most active fields 
in the last decade. The understanding of a number of new and intriguing results, such as the 
biquadratic exchange term in the free magnetic energy of the system, became an exciting 
challenge from both theoretical and experimental point of view. Until recently, it was found 
that the biquadratic exchange coupling was too small when compared to the bilinear term, 
but recent works have proved that it can play a remarkable role in the properties of magnetic 
multilayersEHi. 

It is known that the magnetic properties can depend strongly on the stacking pattern of 
the layers. In this respect, the physical properties of a new class of artificial material, the so- 
called quasiperiodic structures, became recently an attractive field of researc hi. Quasiperi- 
odic structures, which can be idealized as the experimental realization of a one-dimensional 
quasicrystal, are composed from the superposition of two (or more) building blocks that 
are arranged in a desired manner. They can be defined as an intermediate state between 
an ordered system (a periodic crystal) and a disordered one (an amorphous solid), despite 
the purely deterministic rules used to generate them&t They also share a general property, 
which is perhaps their most characteristic one, namely a complex fractal spectra of energy, 
which can be considered as their basic signaturelHil. Their first experimental realization 
(in quasiperiodic GaAs-AlAs heterostructures) was carry out by Merlin and collaborators^, 
and since then they have become a rapidly expanding object of theoretical and experimental 
research. 

Recently, Tsallis and collaborators^ presented a model, based on the most well-known 
and simple deterministic fractal geometry (the triadic Cantor set), to study the thermo- 
dynamic properties. This set is obtained through the repetition of a simple rule: divide 
a given segment into three equal parts, and then eliminate the central one. They showed 
that the specific heat of such a system exhibits a very particular behavior: it oscillates 
log-periodically around a mean value that equals the fractal dimension of the spectrum. 
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Connections with natural spectra showing multifractal properties were also demonstrated 
afterwardsllllil. Furthermore, in a recent publication^, we extended the model described by 
Tsallis et al, to study the specific heat, of real spin waves modes that propagate in quasiperi- 
odic structures, with new features of the specific heat behavior. Throughout these papers, 
the classical Maxwell-Boltzmann statistics were used. 

It is the aim of this work to investigate these properties even further, by studying the 
influence of the biquadratic exchange term in the specific heat spectra of spin waves in 
a magnetic quasiperiodic structure obeying the Fibonacci sequence. Spin-wave studies of 
magnetic systems, where there are contributions of bilinear and biquadratic exchange inter- 
actions, are quite recent&B. Such investigations renewed the interest in spin-wave analysis, 
particularly for non-collinear configurations, since it can give important information about 
the coupling parameters^. The biquadratic coupling in real systems usually favors perpen- 
dicular alignment of the film magnetization, whereas the bilinear one favors both parallel 
or antiparallel configurations. The inclusion of both exchange terms leads to interesting 
physical properties^. 

The plan of this work is as follows. In Section II we present our general theoretical model, 
based on ferromagnetic films described by the Heisenberg Hamiltonian with a biquadratic 
exchange term, which are stacked following a Fibonacci sequence. The multifractal spectra 
then obtained will be used to determine the specific heat spectra for the Fibonacci quasiperi- 
odic magnetic arrangement described in Section III. The numerical results and the discussion 
of their main features are presented in section IV. 

II. THE MODEL 

In this section we follow the spin-wave model of Ref. 20 for quasiperiodic magnetic super- 
lattices with biquadratic exchange coupling. Specifically, we consider magnetic superlattices 
composed of tia layers of material A and ub layers of material B. Materials A and B are 
simple cubic Heisenberg ferromagnets with bulk bilinear exchange constants J a and Jb, and 
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lattice constant a. The spin quantum numbers of the magnetic moments in each material 
are S A and Sb, respectively. In our model we consider that within the bulk of materials A 
and B the magnetic moments interact only by the bilinear exchange couplings J a and Jb- 
However, across the interfaces, the magnetic moments interact by interfacial bilinear (J^) 
and biquadratic (Jbq) exchange couplings. The Hamiltonian for the bulk of each component 
is 



where the sum in the first term is over nearest neighbors j, H is the applied magnetic field 
in the z-direction, and a is equal to A or B. Also, g is the usual Lande factor and hb is 
the Bohr magneton. The spins of sites that are at the interfaces have a Hamiltonian which 
includes an additional biquadratic exchange coupling term, namely: 



It should be remarked that this term is responsible for the new effects found in the specific 
heat reported in this paper. 

The spin wave dispersion relation in a magnetic superlattice is found by considering the 
wave solution inside each material and applying appropriate boundary conditions at the 
interfaces. The solutions for each material can be written as a linear combination of the 
positive- and negative-going solutions of the bulk case, i.e., 



H a = (-1/2)2 J « Si ■ Sj ~ 9fi B H 



(1) 
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These solutions are linked together using the equation of motion for the operator = 
Sf + iSf , which for a site % at an interface and after using the RPA approximation, is (see 
Ref. 20 for details) 

h i s t = g^HoS? + ]T USjSt - SiSf) 

n.n. 

+ E J bA [Sf&Sj - 1) - - 

n.n. 

-[S](2S i - 1) - SjiSt - l)]St }. (5) 



The boundary conditions, after a tedious but straightforward calculation, can be written in 
a matrix form like 



(6) 
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The explicit form of the above matrices can be found elsewhereS Therefore, the explicit 
relation between the Ith and (I + l)th unit cell amplitudes is 



A 



i+i 



v 



l+l 



N^MbN^Ma 



A, 



A' 



(8) 



Here T = N^MbN^Ma is a transfer matrix. This equation, combined with the Bloch 
ansatz yields 



[T-exj>(iQD)\ 



A, 



4 



o. 



(9) 
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The analogous relation between (Ai-i, A'^) and (Ai, A\) combined with equation (9) yields, 

cos(QD) = (1/2)2Y [T] . (10) 

Here Q is the Bloch wavevector of the collective excitation, and D is the size of the super- 
lattice unit cell. Equation (10) follows from the fact that T is an unimodular 2x2 matrix, 
and it describes the bulk modes of spin waves in a magnetic superlattice. Once we know the 
form of the transfer matrix T, the bulk spin wave spectrum is determined. 

Let us briefly review the quasiperiodic sequence considered in this work. First we recall 
the definition of a substitution sequence. Take a finite set £ (here £ = {A, B}) called an 
alphabet and denote by £* the set of all words of finite length that can be written in this 
alphabet. Now let ( be a map from £ to £* by specifying that £ acts on a word by substituting 
each letter (e.g. A) of this word by its corresponding image ((A). A sequence is then called a 
substitution sequence, if it is a fixpoint of (, i.e. if it remains invariant when each letter in the 
sequence is replaced by its image under (. One of the most famous substitution sequences is 
the so-called Fibonacci quasiperiodic sequence, whose substitution rule is A — > ((A) = AB, 
B — > ((B) = A. This sequence can also be constructed by appending the n — 2th generation 
to the n — lth one, i.e., S n = S' n _iS' n _ 2 (n > 2). This algorithm requires the initial conditions 
So = B and Si = A. Some of the first Fibonacci generations are, 

S 2 = [AB] , S 3 = [ABA], S A = [ABAAB] , etc. 

In a given generation S n , the total number of building blocks is given by the Fibonacci 
number F n , which is obtained by the relation F n = F n _i + -F n _2, with Fo = F\ = 1. Also, 
F n _i and -F n _2 are the numbers of building blocks A and B, respectively. As the generation 
order increases (n » 1), the ratio F n /F n _i approaches r = (1 + y / 5)/2, an irrational 
number which is known as the golden mean. From an experimental perspective, a Fibonacci 
superlattice is grown juxtaposing the two magnetic materials A and B according to the 
above described sequence of letters. It can also be shown that the transfer matrices of 
consecutive Fibonacci generations are related by a convenient recursive relation00. 
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T Sn =T Sn _ 2 -T Sn _^ n>2. (11) 

It is easy to see that, from the knowledge of the transfer matrices T So and T Sl , we can 
determine the transfer matrix for any generation and consequently the spin-wave dispersion 
relation. It should be remarked that the matrix Ts 2 recovers the periodic case. 

III. SPECIFIC HEAT SPECTRA 

The spin wave fractal spectra for the Fibonacci superlattices with biquadratic exchange 
coupling is shown in Fig. 1 for a fixed value of the in-plane dimensionless wavevector, namely 
k x a = 2.0. We can clearly see the forbidden and allowed energies versus the Fibonacci 
generation number n, up to their 10th generation, whose unit cell is composed of F 9 = 55 
A's and F 8 = 34 _B's building blocks. The number of allowed bands is equal to three times 
the Fibonacci number F n of the correspondent generation. Notice that, as expected, for 
large n the allowed band regions get narrower and narrower and they have a typical Cantor 
set structure. 

We address now the specific heat of the spectra shown in Fig. 1. The description below, 
which follows the lines of Ref. 17, is general and has been successfully applied to many other 
banded spectra (see also Refs. 23 and 24). In Fig. 1 each spectrum, for a fixed generation 
number n, has m allowed continuous bands. We consider the level density within each band 
to be constant. The partition function for the nth generation is then given by: 

Z n = / p(e)e-*de, (12) 
Jo 

Here (3 = 1/T (by choosing the Boltzmann's constant ks = 1), and we take the density of 
states p(e) = 1. After a straightforward calculation we can write Z n as 

1 2m- 1 

Z n = - R E e-^[l-e-^]. (13) 

P i=l,3,... 

Here the subscript n is the generation number, m is the number of allowed bands and 
= e i+ i — €i is the difference between the top and bottom energy levels of each band. 
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The specific heat is then given by 



which can be written as 



with 



C„(T) = ^[T^], (14) 



CJT) = 1 + ^ - f|, (15) 



2m-l 

k= E k 2 ^-4i^ +1 ]- (16) 

i=l,3,... 



and 



2m- 1 

Sn= E fee^-e^e"^ 1 ]. (17) 

i=l,3,... 

Therefore, once we know the energy spectra of the spin waves which propagates in a 
given sequence's generation of a quasiperiodic structure, we can determine the associated 
specific heat by using (15). 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section we present numerical results obtained for the specific heat of Fibonacci 
magnetic multilayers. In our calculations we have assumed the spin quantum numbers 
Sa = 1-0 and Sb = 1-5, and have taken values of the ratio between the interlayer biquadratic 
and bilinear exchange terms (R = Jbq/Ju) i n accordance with experimental data recently 
measured for different regions of the biquadratic and bilinear exchange couplings' regime0@0. 

Fig. 2 shows the log-log plot of the spin wave specific heat as a function of temperature 
in the low temperature regime and in the absence of the biquadratic term (R = 0). We 
have taken the dimensionless common in-plane wavevector to be k x a = 2.0. As we can 
see, there is an interesting harmonic oscillation of the specific heat, which is a log-periodic 
function of the temperature, i.e., C n (T) = aC n (aT), where a is a constant, and a an 
arbitrary number. These log-periodic oscillations can be traced back to the log-periodicity 
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of the density of state's spectral staircase. They resemble very much the shape of the devil's 
staircase obtained from idealized Cantor sets0'0. The number of oscillations depends on 
the generation number of the Fibonacci sequence: the bigger the system, the greater the 
number of oscillations. Most important, however, is the well defined even and odd parity 
spectra related to the generation number of the Fibonacci structure, with the amplitudes 
of the latter being bigger than the amplitudes of the former. They are due to the parity 
of the set of eigenvalues (e$) used to calculate the partition function given by (13). These 
harmonic oscillations can be identified as the signature of the Fibonacci structure, and it 
has no counterpart in the idealized triadic Cantor set. 

There is an important modification, however, when the biquadratic coupling is present 
(R ^ 0). By contrast with the spectra found in the previous workEl (but with other similar- 
ities remaining) , the two different symmetrical profiles of oscillations are broken for different 
ratios between the biquadratic and bilinear terms. To reinforce this intriguing behavior, 
we show, in Figs. 3 and 4, the evolution of this broken-symmetry for R = 0.2 (where the 
even-mode symmetry starts to be broken ) and 0.4 (where the odd-mode symmetry also 
starts to be broken), respectively. A possible explanation is that these different behaviors 
arise because the biquadratic exchange term reduces the effective coupling between the ad- 
jacent magnetic layers at the interface, e.g., as may be deduced by examining the overall 
coefficients of or Sj~ in (5) (see also Ref. 20). Also the biquadratic exchange induces 
long range correlations that emphasize the quasiperiodicity of the system. In other words, 
it looks like these effects make the whole structure see better its quasiperiodicity, in effect 
increasing its degree of disorder as R increases! It turns out that the spin wave energy 
band structure is more disordered and, as a consequence, this is reflected in the specific heat 
curves, which do not have the well-defined standard profiles (actually the oscillations are 
non-harmonic!) found in the absence of the biquadratic term, as well as presented for other 
excitations^'^. This argument is reinforced by previous works on the correlation lengths of 
magnetic systems exhibiting biquadratic exchange coupling (see, for example, S0rensen and 
Youngil). 
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To summarize, we have proposed in this paper a realistic model to study the specific heat 
contribution from the energy spectra of real spin waves in quasicrystals of the Fibonacci type, 
stressing the important role played by the biquadratic exchange term. Certainly the theoret- 
ical predictions shown here can be tested experimentally, and we encourage experimentalists 
to carry out such studies. 

Further investigations for the case of the specific heat due to propagation of mixed modes 
is currently under consideration, and we hope to present these results in a later publication. 
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Figure Captions 



1. Energy spectra of spin waves for the quasiperiodic Fibonacci structure when a bi- 
quadratic interlayer exchange term is present. Here we consider the dimensionless 
in-plane wavevector k x a = 2.0. 

2. Log- log plot of the specific heat versus temperature for the generation numbers of the 
Fibonacci quasiperiodic sequence in the absence of the biquadractic term: (a) even 
generation numbers; (b) odd generation numbers. 

3. Same as Fig. 2, but for the ratio between the interlayer biquadratic and bilinear ex- 
change terms R = 0.2. From this value the even symmetry starts to be broken. 

4. Same as Fig. 2, but for the ratio between the interlayer biquadratic and bilinear ex- 
change terms R = 0.4. From this value the odd symmetry also starts to be broken. 
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